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Abstract 

Dirac showed that the existence of one magnetic pole in the universe could offer 
an explanation of the discrete nature of the electric charge. Magnetic poles appear 
naturally in most grand unified theories. Their discovery would be of greatest impor- 
tance for particle physics and cosmology. The intense experimental search carried 
thus far has not met with success. I proposed a universe with magnetic poles which 
are not observed free because they hide in deeply bound monopole-anti-monopole 
states named monopolium. I discuss the realization of this proposal and its consis- 
tency with known cosmological features. I furthermore analyze its implications and 
the experimental signatures that confirm the scenario. 
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1 Introduction 



The theoretical justification for the existence of classical magnetic poles, hereafter 
called monopoles, is that they add symmetry to Maxwell's equations and explain 
charge quantization [Tj . Dirac showed that the mere existence of a monopole in the 
universe could offer an explanation of the discrete nature of the electric charge. His 
analysis leads to the so called Dirac Quantization Condition (DQC), 

^ = - N = i2 m 

where e is the electron charge and g the monopole charge [1]. Note that if quarks 
were asymptotic states the minimum monopole charge would be three times larger. 

The origin of monopoles, and therefore their properties, is diverse. In Dirac's 
formulation monopoles are assumed to exist as point-like particles and quantum 
mechanical consistency conditions lead to Eq.(l), establishing the value of their 
magnetic charge. However, their mass, m, is a parameter of the theory, limited 
only by classical reasonings to be m > 2 GeV [2]. In non-Abelian gauge theories 
monopoles arise as topologically stable solutions through spontaneous breaking via 
the Kibble mechanism [3]. They are allowed by most Grand Unified Theory (GUT) 
models, have finite size and come out extremely massive m > 10^^ GeV. Further- 
more, there are also models based on other mechanisms with masses between those 
two extremes [2], HI |5] . 

The discovery of monopoles would be of greatest importance not only for particle 
physics but for cosmology as well. Therefore monopoles and their experimental de- 
tection have been a subject of much study since many believe in Dirac's statement pQ 

"...one would be surprised if Nature had made no use of it [the monopole] 

At present, despite intense experimental search, there is no evidence of their 
existence [21 HI EJ [TJ [S]. This state of affairs has led me to investigate a possible 
mechanism by which monopoles could exist and still be undetectable by present 
experiments. 

Although monopoles symmetrize Maxwell's equations in form there is a numerical 
asymmetry arising from the DQC, namely that the basic magnetic charge is much 
larger than the smallest electric charge. This led Dirac himself in his 1931 paper [T] 
to state, 

"... the attractive force between two one quantum poles of opposite sign is 
^i|7-j2 ^ 4692| time that between the electron and the proton. This very large 
force may perhaps account for why the monopoles have never been separated." 

This statement by Dirac has motivated the present investigation. I propose a 
scenario where monopoles exist but are hidden from our direct observation because 
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today they appear forming deep monopole-anti-monopole bound states. I introduce 
in the next paragraph my proposal which I elaborate in detail in the next sections 
together with its experimental connotations. 

At some early stage in the expansion of the Universe, monopoles and their an- 
tiparticles were created. At a later time the dynamics was such that most of the 
poles paired up to form monopole-anti-monopole bound states called monopolia. 
This happened because it is easier for the monopoles to interact with each other 
than with the charged particles in the hot plasma, thus as a consequence also no 
significant friction force arises. Therefore, the lifetime of the primordial monop- 
olium is solely governed by cascading to the lower bound states where the poles 
finally annihilate. Moreover, the lifetime of monopolium is sufficiently long to allow 
primordial monopolia to exist even today in measurable abundances. Thus today, 
almost all of the existing monopoles, appear confined in deeply bound states [9]. 
However, the present investigation would be irrelevant if no proof of the existence 
of monopolium could be found. I analyze signatures of their present existence and 
of their formation period. 



2 Hidden monopoles 

I envisage a scenario, in which monopoles are not observable as free states at present, 
which is realized by means of a few assumptions, that satisfies all phenomenological 
restrictions and leads to new observations which can sustain it. If my scenario is 
confirmed experimentally it will strongly restrict the way cosmological models deal 
with monopoles. 

At some early stage in the expansion of the universe monopoles and their an- 
tiparticles were created by a mechanism which is free from the standard monopole 
problem [TOj. No precise mechanism for their creation is advocated, therefore the 
mass is not fixed and is left as a parameter to be fitted by consistency requirements. 
Moreover, monopoles and anti-monopoles existed in the universe, at that time, at 
the level of abundance compatible with known phenomenological and experimental 
upper bounds [HI [T2|, [13]. These are the same starting assumptions of all similar 
treatments [HI [151 [IS]. I next depart from them by assuming that during nucle- 
osynthesis, most of the monopoles and anti-monopoles bind in pairs, due to the 
strong magnetic forces, to form monopolium. This scheme is realized physically by 
imposing that: i) the capture radius of the poles, r capture, and the mean free path 
for charge particle collision in the hot plasma. A, satisfy, 

^capture A, (2) 

ii) a consistent monopolium formation scenario. Let me discuss in this section the 
first assumption which leads to the determination of some of the monopole and 
monopolium properties and leave for the next section the discussion on monopolium 
formation. 
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The capture radius, r capture, is given by [21 [9], 



f capture ^ j^rp • (^) 

The mean free path A is given by 

A ~ ^ , (4) 
where a is the cross section for the scattering of the monopole with charged particles 

mm 

2 

TTIC 

CT ~ 2 — — nanobarns , (5) 

Pch is the density of charged particles and m is the mass of the monopole. 

I describe the monopoles from the monopolium formation era up to the present 
days by point like Dirac monopoles. It is reasonable to do so since the discussion is 
largely independent of the detailed structure of the monopoles because it depends 
only on global properties, i.e., magnetic charge, mass and cosmological abundances. 

In the rest of the paper I proceed to show that my assumptions lead to a picture 
which is consistent with present data. 

I describe monopolium as a Bohr atom, with reduced mass m/2 and a strong 
magnetic, instead of a weak electric, coupling. Its binding energy is 

\8aJ 

where a = ^ i?> the fine structure constant of QED and n the principal quantum 
number. This equation and those that follow are to be considered only for large 
principal quantum number {n > 50). For low values of n the annihilation mechanism 
becomes dominant. 

The approximate size of the system is given by 

2 /^X 

r ~ < r >„ ~ an. u) 

mc 

To calculate the mean life I distinguish two processes, i) the cascading process, 
dominated by dipole radiation which I apply as 

T,^poie - ~ 2 (12)=^ ^ a' nt , (8) 

where Ui is the principal quantum number associated with the initial bound state 
which will be very large Ui ~ 10^; ii) the annihilation process, which due to the 
magnitude of g is highly non perturbative and which I next estimate. Looking at 
the two photon decay process I see [TTl US] 
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Tannihilation << Tg^ ~ 2 (4)^ tt^ , (9) 

where Uf is the largest principal quantum number associated with a state at which 
annihilation is still efficient. Since the monopole and anti-monopole only annihilate 
efficiently when there is a considerable probability of being on top of each other and 
this only happens for n < 50, Uf << rii. Thus the annihilation mean life is small 
compared with the cascading time and can be disregarded in the time scale analysis. 

The previous equations can be summarized in terms of the binding energy of the 
initial bound state and the mean life as, 



Eb{eV)n{A) ~ 5.10^eVA (10) 
~ 9. 10^[Eb(e\/)i/V(sec)i/^] (11) 
mc^eV) ~ 3. 10'[Eb{eVf/^T{secy/^] eV . (12) 

Here Eb and r?, are respectively the binding energy and the radius of the initial 
bound state. 

Using Eqs. ([2]) through f|T2l) my scenario can be constructed. I assume that cap- 
ture takes place for a binding energy slightly higher, to avoid thermal dissociation, 
than kT = 1 MeV, 

Eb > 1 MeV . (13) 

This temperature is not related to the scale for production of monopoles but to that 
at which monopolium, the bound state, is formed from already existing monopoles 

m- 

These equations and the temperature scale of the proposed scenario show that mo- 
nopolium is a very tightly bound system 

n < 0.05A < r^apture ~ O.OTA . (14) 

Let me now proceed to the calculation of the mean free path, 

A ~ . 

(^Pch 

The cross section at kT = 1 MeV is given by 

cr ~ 2. 10"^^ ^ — - A . 

eV 

From Eq.(12) I have 

mc^feV) 1 

^—^ = 1.5 10^^ r/2 , 

eV 
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where rj = t/tu, t being the hfetime of the monopohum and tjj the age of the 
universe, which we have taken approximately to he tu = 1/Hq 3. 10^'' sec. 
Thus the cross section becomes 



(T ~ 3. 10^ ri^ 



(15) 



I now calculate the density of charged particles. The density of photons is given 



while the ratio of the densities of nucleons to photons is given in nucleosynthesis 
by Pn/Pi ~ 4. — 7. 10~^°. At kT = 1 MeV , we have mostly nucleons and since 
pnl Pp ~ exp (— 1.293//cT), we obtain 



In Fig. [T]I show the ratio of A over r capture as a function of the life of monopolium 
in units of the age of the universe. I also plot the ratio of the average distance 
between charges to Vcapture- The plot shows that for lifetimes above 0.1 tjj the mean 
free path is smaller than the mean distance between charges and the capture radius 
becomes of the order of magnitude of the mean free path. Thus for monopolium 
lifetimes of the order of the age of the universe the conventional scenarios would take 
place [m [Tbilisi- However, for t/tjj < 0.1 the condition Eq.(2) will be satisfied and 
a completely different scenario occurs whose existence has been a matter of debate 
since the idea was first proposed [19] . 

In order to clarify the issue let me recall once more the work of Blanco-Pillado 
and Olum [TB]. They pointed out that in the standard scenario the drag force felt 
by the monopole in the plasma reduces dramatically the mean life of the state [TS| 
and therefore relic monopolia do not influence present day observations. This is so 
if we are to the left of t/tu > 0.1. However, their phenomenon is not active in the 
scenario presented here, i.e., when t/tu < 0.1 for obvious reasons. As shown in 
ref.[2T] the drag force depends on the two limits of an impact parameter integral. 



by m 



2C(3)(A:T)3 




(17) 



(16) 
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Figure 1: The ratio \/r capture is presented as a function of monopolium lifetime r in units 
of age of the universe tu at kT = 1 MeV. Also the ratio of the average distance between 
charge particles < dch > /f capture at the same temperature is shown. 

where naturally the force only arises if 6mm < ^max- The small limit has to do with 
the monopole interaction with the plasma, in our case hmin ~ A. The large limit 
depends on the interaction of the monopoles inside the monopolium with the charge 
particles and can not extend beyond r captures bmax ~ f capture- Due to our assumption 
Eq.(I2]), as seen in Fig. [H r capture < A in the region of interest and therefore in 
our scenario there is no drag force. It is interesting to note that the solution that 
Blanco-Pillado and Olum find to solve the impasse, i.e. monopoles attached by 
strings [16], reproduce with a complex dynamics the same scenario for monopolium 
I obtain, i.e. non-relativistic monopoles very closely tight together so that the light 
charges do not interact with them and they only loose energy by radiation. 

3 Monopolium formation 

In the previous section I have shown that the proposed scenario might be realized 
if ?7 << 1 since under these conditions primordial monopolium might be observed 
today. I here study the other ingredient of the scheme, a plausible scenario for the 
formation rate of monopolium which gives consistency to the scheme. The scheme is 
based on three time (temperature) scales. The first time scale is the formation scale 
of monopoles, tj, whose corresponding temperature is Tj. This temperature is very 
high, consistent with the large mass of the monopoles determined in the previous 
section. The second time scale ends at the beginning of the monopolium bound 
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state stability period, i.e. when the temperature is lower then the minimum binding 
energy, i.e. Tf ~ IMeV. Finally the third time scale is today, i.e. Ttoday ~ 2.7 ^K. 

In this section I will study the formation rate of monopolium from the early 
universe to the beginning of the stability period, i.e. from Tj to Tf. In the next 
section I will study the decay process to find present day monopolium abundance, 
i.e. from Tj to Ttoday 

Following Blanco-Pillado and Olum [16] one can write the evolution equation as 



r/r / 1 \ (9/^°) 



In this equation F represents the comoving monopolium density, i.e. N]^jfj/ s, where 
^MM is the monopolium density and s the entropy density; 7 is the monopole 
comoving density, i.e. n^/s, where um is the monopole density; T is the temperature 
and A an softly temperature dependent quantity. 
Let 7 depend on temperature as 

7 = BT\ (19) 

where B is a constant, a possibility contemplated in ref. ^0]. The above equation 
can be easily solved leading to 

F(T.) - T{Tj) = -^^^ (T.) Tt' - i\Tj) Tf') (20) 
My fundamental hypothesis implies. 



F(T,) « F(T;), 
7(T,) » 7(T^). 



Therefore, 

A 



T{Tf)^,^^^^\T,)Tf\ (21) 
Dividing by 7(T/) I arrive at, 

NMMiTf) ^ AB (Tl\ ' 

UMiTf) ~ 25 + 0.1'^ [tj] ■ ^''^ 

Thus if 6 > 0, since Tj >> Tf, the wishful scenario is realized. Note that this 
behavior corresponds to p > 1 in ref. j20j . 
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4 Monopolium abundance 



I proceed in here to investigate phenomenological consistencies of the proposed pic- 
ture. Let us calculate the present abundance of monopolium taking as input the 
abundance during the formation period. By doing so my aim is to find consistency 
between the proposed scenario and the observation. 

The equation governing the density of monopolia p taking into account the ex- 
pansion of the universe is given by 

p{t) = -\ Pit) - 3 I Pit) , (23) 

where R is the scale factor of the universe and A denotes dA/dt. If the expansion is 
adiabatic (i?T ~ constant) and the universe is radiation dominated, the expansion 
rate is given by 

^ = = (24) 

R T A ^ ' 

Here, T is the temperature of the universe and A a quantity related to the Planck 
mass and the effective degrees of freedom [20j. The temperature equation in Eq. fl2^ 
can be integrated to give 

t = ^ (25) 

establishing the relation between evolution time and temperature. Using Eqs. fl2^ . 
(124|) and ([25]) one obtains, 

Pit) = -(z + ll)pit) (26) 



r 2 1 



which can be easily integrated giving 



Pit) = Pito) {--^) ■ (27) 

I now take this equation, which is the conventional equation for the decay of un 
unstable system in an expanding universe pH], and adapt it to my interests, namely 
I want to study the evolution of the number of monopolia between Tf and Ttoday 
Using that, in this case, t = ttoday = W » to = t(T/) ~ 1 sec, and the relation 
between time and temperature, Eq. (125|1 . one can write 

piTtoday) = piTj) (^y^'exp (-^) . (28) 
In the present situation this equation reduces to 

p(2.7°ir) = 1.25 10-2^ pikT = IMeV) exp i--) , (29) 

V 
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Figure 2: The density in a cm^ {logio) of monopolia (solid curve) and the number (logio) 
of monopoha decays in a pc^ (dashed) today presented as a function of monopohum hfetime 
T in units of age of the Universe tjj . The values by Hill \2A\ (dot) and Blanco-Pillado 
and Olum |16j (short dot) for the density of monopolia are used to extract the number of 
decays. 



where rj = t/tu. 

The standard scenario for helium synthesis requires that the mass of the monopole 
does not dominate the universe when kT = 1 MeV, this implies [20] , 

mc^ (eV)\ 



PmonopoieikT = IMeV) < 1.3 10^^ ^ A-^^ 



I assume that the density of monopoles is equal to the density of anti-monopoles 
and since most of them are bound equal to the density of monopolia, thus one gets 
from Eq. fl2^ for the density of monopolia today, 

p(2.7°K)^1.6 10-^ (a^lm]" A-s^ 



eV 



which using Eq. (fT2l) leads to 



p(2.7°K) ~ lO-^*^ """"^[y cm-^ . (30) 

This equation establishes a relation between the parameter characterizing the various 
scenarios rj and the observation. 

Let me calculate the number of decays per year in a given volume of the universe 
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Figure 3: The density in a cm^ (^0510) of monopolia (solid curve) and the number in a 
pc^ (logio) of monopoha decays (dashed) today presented as a function of the monopolium 
hfetime r in units of age of the Universe tjj . The hmits of the values for the density of 
monopolia from the values of Bhattarchee and Sigl [15], 1 < ^Mh'^ < 10^ are plotted and 
the corresponding number of decays shown. 



. = iv(l-exp(-i))~^. (31) 

where I have used the approximation that riyears) » 1 year. Note also that 
the temperature factor in Eg. (1231) drops out because it is very close to 1. is the 
number of monopolia in the given volume. Let me choose to calculate the observation 
a volume of 1 (pc)^, then I get for the number of events in one year 

u ~ 2.7 10^9 ^^^^^1^ . (32) 

Monopolium has been associated with Ultra High Energy Cosmic Rays (UHECR) [221 
123] in various schemes [HI [TSl US] ■ This association leads to a phenomenological de- 
termination of its abundance. I look for consistency between the phenomenological 
determined abundances and the monopolium mean life obtained in the calculation. 

In Fig. (2) we show the results of our calculation as a function of r in units of 
Tjj and I introduce two numerical values for the density as obtained by 

i) Hill [H]: p(2.7°K) ~ IQ-^^cm^^ 

ii) Blanco-Pillado and Olum [16] : p(2.7°K) ~ lO'^^cm-^ 
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In Fig. (3) I compare my results to bands, characterized by 1 < VLm^? < 10^, 
obtained from the density equation of Bhattacharjee and Sigl [15]. 

The above study shows that proposed values for the density of monopolia are 
consistent with t] ~ 0.02 and therefore with 

A ~ 10 r capture 

If we look at monopolia decays we realize that their number is 

10^ < z/ < 10^^ decays/pc^ 

which is large compared with those obtained in other models. 

Please note that all the estimates are based on the assumption that all UHECRs 
are due to monopolia. If they are additional mechanisms for UHECRs formation 
or if the number of observations [22l [23] diminish, the density of monopolia today 
would decrease and the proposed scenario would become more natural. However, 
this would also imply that it would be harder to confirm experimentally. 

I conclude from the above analysis that at present most monopoles in the inter- 
galactic vacuum are bound in deep bound states n ~ 50 close to the annihilation 
levels and therefore their binding energy is at the level of 

~ 10^^ GeV , (33) 

supporting Dirac's conjecture for the non observability of monopoles. 
Moreover, by looking back at Eq.(12) one sees that 

mc^ ~ 2. 10^^ GeV, (34) 
which is very large and comparable to the values arising from GUT models. 

5 Monopolium detection 

Ideally we would like to be able to produce the monopolium in our laboratory. Its 
mass M ~ 10^^ GeV makes laboratory production impossible. Could we capture 
a monopolium in our laboratory and measure its properties? It soon will become 
clear this is an impossible task. Monopolium is a sterile particle under laboratory 
conditions. 

The present day background monopolium density is small. 

p(2.7°K) < 10"^^ cm-l 
At present monopolia are mostly thermal and therefore their velocity is of the order 

V ~ 10~^ m/ sec. 
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Thus the average number of particles in our detector would be 

Nm ^ pvt A < 10"^^ (t/years) (A/Km^), 

where A is the area of the detector and t the time of exposure. Thus the chances of 
having one localized event are insignificant. 

Let me assume that a monopolium enters by chance a detector, what could we 
observe? Due to the dual behavior of the Maxwell equations in the presence of 
monopoles [21], monopolium has an electric dipole moment 

9 f 

P ~ 7^ L 

2 m c 

and a magnetic dipole moment 

flr^2gc< r>i^o 

In the presence of an electric or magnetic field they tend to orient against the field 
and become effectively 

9 

P ~ m 

2 m c 

and 

c 

/i~ 2 --(3^2 -/(/ + !)), 
m g 

where n, I, mi are respectively the principal, the orbital and the magnetic quantum 
numbers. Substituting the values of the couplings and masses I get 

p ~ lO'^^m/ (e-charge) meter 

and 

/i ~ 10~^2 (3^2 - / (/ + 1)) eV/Tesla 

Note that, the maximum achievable electric and magnetic field gradients today , 50 
MV/meter^ and 5 Tesla/meter, produce in background monopolia, with n ~ Z ~ 
mi ~ 100, insignificantly small energies and forces. It is easy to calculate, that 
to stop them, moving at thermal velocities, we would need distances and times of 
universe scales. 

Thus I conclude, that from the point of view of traditional laboratory experi- 
ments, monopolia are sterile and we have to center our attention in astrophysical 
observations. 

Let me now turn to astrophysical observations. Let me distinguish two periods: 
the formation period and the immediate past. 

During the formation period n ~ 10^. Cascading for large values of n leads to 
Larmor type emission 
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A ~ 16a^—nl ~ A , (35) 

mc mc^ 

thus the wave length during the nucleosynthesis period will be 

A ~ milimeter. 

Therefore, there should be an isotropic background radio frequencies as a remnant 
of that period. 

In the immediate past even until today, most monopolia are close to the annihi- 
lating stage, i.e., n < 100. The Larmor formula is still approximately applicable, 
thus 

A < 10"^^ A, 

which implies that the emitted photons will have a huge energy. 

hu ~ 10^° GeV 

This energetic photons will occur at the level of 10^^ — 10^° per pc^ and therefore 
should be seen. However, their distribution will not be isotropic. The core of galax- 
ies, and of clusters of galaxies, provide an environment of high electric and magnetic 
field gradients, thus the small electric and magnetic dipoles of monopolium will 
change their distribution in these environments and create geographic anisotropies. 

Moreover, these regions also provide an environment with high energy and high 
density where monopolia might be excited to the point of break up. Thus, also 
the low frequency spectrum will acquire geographical distributions which become 
anisotropies in the spectrum over the isotropic spectrum remnant of the formation 
era. Note that very few monopolia can be formed after electron-positron annihilation 
due to the lack of remnant monopoles and the absence of monopole pair formation 
due to their huge masses, except in these very energetic environments, and this 
phenomenon will not affect greatly the calculated density of monopolia [l6j. 

Finally, monopolium can annihilate into UHECRs, at the level of millions per year 
and per cubic parsec, under present experimental expectations [221 123] provided it is 
the only mechanism, depositing a huge amount of energy, E > 10^^ GeV, in a small 
region of space-time leading to what Hill [T3] calls a cataclysmic scenario, whose 
details depend on the microscopic theory of monopole formation. 



6 Conclusions 

The possibility of having monopoles in nature is appealing. 1 have presented a 
scenario for the universe in which relic monopoles still exist today however, not as 
free particles, but deeply bound in monopolium states. The crucial ingredients of my 
proposal are: i) that, in the early universe, the mean free path of monopoles is much 
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larger than their capture radius, and therefore they bind so tightly in monopohum 
that they barely interact with the surrounding plasma, surviving in this way the 
effect of the drag force and only emitting energy by radiation until their annihilation; 
ii) the evolution of the monopole density is governed by a strongly temperature 
dependent function leading to a large production of monopolia. Few monopolia are 
formed in the second period since almost no free monopoles exist, because they are 
to produce at these low temperatures due to their large mass and strong binding 
within monopolia, to drive the formation rate. The initial density of monopolia and 
their lifetime might explain UHECRs by construction. 

Three distinctive quantities determine the consistency of the various requirements 
in the scheme: the monopole evolution parameter 6; the temperature of monopo- 
lium formation {kT ~ 1 MeV); and the mean life of the state (r ~ 10^ years). 
The outcome is monopolium, a neutral particle protected from the interaction with 
the medium in a strongly bound state, which radiates copiously until ultimately 
annihilates in a cataclysmic scenario being a possible source UHECRs. 

The detection of monopolia, and therefore the existence of monopoles, presents 
interesting signatures associated with the monopolium spectrum, i.e., a diffuse 
isotropic radio frequency background remnant of its formation period with geo- 
graphical anisotropics produced by more recent activities, and a high frequency mo- 
nopolium spectrum associated with the last period of its lifetime, which manifests 
itself by anisotropically distributed high energy gamma rays. 
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